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I.  RESULTS  FROM  PRIOR  SUPPORT 

1.  Simulation  Study  of  a  New  Mechanism  for  Excitation  of  Kinetic  Waves  in  a 
Magnetoplasma 

Nishikawa  et  al.  [1988]  have  investigated  the  new  ion-cyclotron-like  waves  by  a  localized 
transverse  electric  field  by  means  of  simulation  with  the  assistance  of  the  nonlocal  kinetic 
theory.  The  linear  theory  shows  that  the  growth  rates  of  the  kinetic  Kelvin-Helmholtz  (K-H) 
modes  are  strongly  reduced  with  increasing  u  =  k\\/ky,  and  they  become  unstable  only  when 
b  =  kyp]  <  1  and  kyL  ss  1.  where  L  is  the  scale  length  associated  with  the  transverse  electric 
field.  On  the  other  hand,  the  new  modes  have  larger  frequencies  and  become  unstable  at  larger 
b  >  1  and  kyL  >  1  [Ganguli  et  al.,  1988]. 

Simulation  results  show  that  ion-cyclotron-like  waves  are  excited  in  regions  where  E  x  B 
drift  is  localized.  Linear  growth  rates  of  several  modes  are  obtained  from  the  wave  analysis 
of  the  simulation.  This  linear  analysis  shows  that  the  (0,  4)  mode  corresponds  to  large  b.  and 
large  kyL  has  the  maximum  growth  rate.  Clearly,  these  are  not  K-H  modes.  Further,  the 
simulation  results  show  that  density  gradients  help  to  enhance  the  growth  rates.  However,  like 
the  K-H  mode,  the  real  frequencies  of  this  instability  are  approximately  proportional  to  ky Vg. 
where  Vg  is  the  peak  value  of  the  E  x  B  drift. 

Nonlinear  phenomena  such  as  diffusion  and  coalescence  of  vortices  are  investigated.  In  the 
linear  stage  smaller  vortices  are  generated  and  larger  vortices  with  the  lower  frequencies  are 
dominant  in  the  nonlinear  stage.  In  the  nonlinear  stage  ions  diffuse  strongly  due  to  large-scale 
vortices. 


1 


Recently,  we  have  investigated  the  electrostatic  waves  driven  by  the  combined  effects  of  a 
localized  transverse  electric  field  and  parallel  electron  drifts  by  means  of  simulation  with  the 


assistance  of  the  nonlocal  kinetic  theory  [Nishikavva  et  al.,  1989b]. 

We  have  performed  a  number  of  simulations  for  this  instability.  Simulation  results  show 
that  electrostatic  waves  are  excited  in  the  regions  where  the  E  x  B  drift  is  localized  in  the 
simultaneous  presence  of  parallel  electron  drifts  and  transverse  electric  fields.  Simulations  with 
only  parallel  electron  drifts  or  transverse  electric  fields  show  that  no  instability  grows  out  of 
the  thermal  noise.  The  simulation  with  both  the  parallel  electron  drift  and  the  transverse 
electric  field  shows  the  growing  waves  out  of  the  thermal  noise.  The  Doppler  shift  due  to  the 
E  X  B  drift  can  lower  the  phase  velocity  of  waves  along  the  magnetic  field.  Then  this  makes 
wave-particle  resonance  possible  for  smaller  V£,  which  leads  to  this  instability. 

The  nonlinear  phenomena  such  as  diffusion  and  coalescence  of  vortices  are  investigated.  In 
the  linear  stage,  smaller  vortices  are  generated  and  larger  vortices  with  the  lower  frequencies 
are  dominant  in  the  nonlinear  stage.  In  the  nonlinear  stage  the  ions  diffuse  strongly  due  to 
the  large  scale  vortices. 
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2.  Beam  Instability  in  the  Foreshock 

As  an  application  of  the  simulation  method  used  in  the  proposed  research  (Broadband 
electrostatic  noise),  the  beam  instability  in  the  foreshock  has  been  investigated.  Electrons 
backstreaming  into  the  Earth’s  foreshock  generate  waves  near  the  plasma  frequency  by  the 
beam  instability.  Two  versions  of  the  beam  instability  exist:  the  'reactive'  version  in  which 
narrowband  waves  grow  by  bunching  the  electrons  in  space,  and  the  ‘kinetic’  version  in  which 
broadband  growth  occurs  by  a  maser  mechanism  [Cairns,  1887a,  b,  and  references  therein].  Re¬ 
cently,  Cairns  [1987b]  has  suggested  that  (1)  the  backstreaming  electrons  have  steep-sided  'cut¬ 
off’  distributions  which  are  initially  unstable  to  the  reactive  instability,  (2)  the  back-reaction 
to  the  wave  growth  causes  the  instability  to  pass  into  its  kinetic  phase,  and  (3)  the  kinetic 
instability  saturates  by  quasilinear  relaxation. 

Cairns  and  Nishikawa  [1989]  have  performed  two-dimensional  simulations  of  the  reactive 
instability  for  Maxwellian  beams  and  cutoff  distributions.  The  results  of  the  simulations  are 
consistent  with  suggestions  (1)  and  (2)  above.  In  addition,  we  have  demonstrated  that  the 
reactive  instability  is  a  bunching  instability,  and  the  reactive  instability  saturates  and  passes 
over  into  the  kinetic  phase  by  particle  trapping.  We  found  that  the  kinetic  growth  occurring 
after  saturation  of  the  reactive  instability  is  presumably  due  to  the  spatially  localized  gradients 
in  y  —  v ||  phase  space.  Both  simulation  results  and  numerical  solutions  of  the  dispersion  equation 
indicate  that  the  center  frequency  of  the  intense  narrowband  waves  near  the  foreshock  boundary- 
may  be  between  0.9cupe  and  0.98wpe,  rather  than  being  above  uipe  as  previously  believed. 
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3.  Whistler  Mode  Driven  by  the  Spacelab  2  Electron  Beam 

During  the  Spacelab  2  mission  while  an  electron  beam  was  being  ejected  from  the  shuttle, 
the  Plasma  Diagnostics  Package  (PDP)  detected  a  clear  funnel-shaped  emission  that  is  believed 
to  be  caused  by  whistler-mode  emission  from  the  electron  beam  [Gurnett  et  ai,  1986].  In  order 
to  understand  the  mechanism  of  this  emission,  the  simulations  have  been  performed  using 
a  three-dimensional  magnetostatic  code  for  low-/?  plasmas  in  which  the  beam  electrons  are 
initially  located  in  the  column  [Nishikawa  et  al. ,  1989a].  In  order  to  simulate  the  continuous 
electron  ejection  from  the  shuttle,  the  simulations  were  also  performed  with  the  recycling  of 
the  beam  electrons.  The  beam  electrons  excite  whistler  waves  and  lower  hybrid  waves.  The 
brief  fluid  theory  based  on  the  magnetostatic  code  was  checked  with  the  simulation  results. 
The  propagating  whistler  mode  was  identified  with  the  theory.  The  simulation  results  show 
that  the  quasi  perpendicularly  (the  angle  between  the  magnetic  field  and  the  wavenumber  is 
larger  than  50°)  propagating  whistler  waves  have  larger  amplitude  whose  real  frequencies  are 
smaller  than  the  local  electron  cyclotron  frequency.  This  fact  is  consistent  with  the  fact  that 
the  funnel-shaped  emission  is  observed  below  the  electron  cyclotron  frequency  away  from  the 
beam  electron.  The  beam  electrons  initially  in  the  column  diffuse  radially  as  well  as  slow  down 
due  to  the  E  X  B  caused  by  the  excited  waves.  The  acceleration  of  the  beam  electrons  also 
takes  place  due  to  the  excited  whistler  waves. 

In  order  to  compare  with  the  PDP  data,  the  local  magnetic  field.-.  Bx,y  and  the  perturbed 
electric  fields  Ey<z  are  diagnosed  at  the  several  points  in  the  simulation  system.  The  results 
show  that  the  waves  are  radially  excited  by  the  beam  electrons  localized  in  the  center  of  the 
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system.  The  wave  spectra  of  the  electric  fields  Ey,z  diagnosed  at  x  =  3 1  Ae,  y  =  16Af.  and 
s  =  17Ae  show  the  several  kinds  of  waves  generated  by  the  beam  electrons.  The  analysis 
shows  that  the  lower  hybrid  waves  and  whistler  waves  are  excited.  The  frequency  range  of 
these  spectra  extends  from  u>pt  to  beyond  u>pe  which  is  in  qualitative  agreement  with  the 
PDP  data.  As  the  PDP  data  show,  the  intense  electrostatic  narrowband  emission  around  the 
electron  plasma  frequency  has  been  observed  dominantly  in  the  spectra  of  the  electric  field  ( Ez ) 
with  and  without  the  recycling  of  the  beam  electrons.  This  wave  is  identified  as  the  parallel 
(quasi  parallel)  whistler  wave.  This  simulation  result  is  in  good  agreement  with  the  PDP  data 
[Gurnett  et  al.,  1986].  The  spectra  of  the  magnetic  fields  diagnosed  at  the  same  position  show 
that  whistler  modes  are  excited.  The  cB/E  was  calculated  for  the  whistler  mode.  In  the  case 
without  the  recycling  of  the  beam  electrons,  the  cB/E  is  approximately  0.35  for  the  whistler 
mode  (around  0.368wpP),  which  is  smaller  than  that  calculated  from  the  PDP  data  [Gurnett 
et  al.,  1986]. 
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